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A QCD fit analysis to tiie combined HERA inclusive deep inelastic cross sections measured by 
the HI and ZEUS collaborations for e p scattering to extract HERAPDF sets is presented. 
The results are used for predictions of pp processes at Tevatron and pp processes at the 
LHC. The QCD analysis has been extended to include the combined HERA 11 measurements 
at high Q^ resulting in the HERAPDF1.5 sets, with full estimation of uncertainties. The 
precision of the new PDFs at high x is considerably improved, particularly in the valence 
sector. In addition, inclusion of the HERA jet data allows for a precise determination of the 
strong coupling. Moreover, inclusion of the preliminary combined HERA charm data provides 
constraints for the optimal value of the charm mass used in QCD theory models which may 
account for some of the differences among global PDF fits. 



1 Combined HI and ZEUS Cross Section Measurements 

The main information on proton structure functions comes from the Deep Inelastic Scattering 
(DIS) collider experiments HI and ZEUS at HERA. Measurements at HERA go well beyond 
the phase space accessible by fixed target experiments with an extended kinematic range of 
0.045 < Q^ < 30 000 GeV^ and 6 x 10^^ < x < 0.65. To further benefit from the precision of the 
measurements the HI and ZEUS cross sections are combined, resulting in the most consistent 
and precise DIS inclusive double differential cross-section measurement of neutral and charged 
current e^p scattering to datel^. The data combination "procedure" has been repeated to 
include as well the preliminary high precision measurements at high Q^ from the second run 
period of HERA. Therefore, these accurate measurements can be used as the sole input to the 
QCD analysis to determine the proton parton distribution functions (PDFs) as described in the 
following sections, which can be used then for precise predictions for LHC processes. 



2 QCD Analysis settings 

HERA PDFs are determined from QCD fits to HERA data alone. Only the region where 
perturbative QCD is valid, data above Q^in = 3. 5 GeV^ are used in the central fit. The HERA 
data have a minimum invariant mass of the hadronic system, W, of 15 GeV, such that they are 
in a kinematic region where there is no sensitivity to non-perturbative effects common to fixed 
target data. 

The fit procedure starts by parametrising PDFs at a starting scale Qq = 1.9 GeV^, chosen 
to be below the charm mass threshold. The QCD settings are optimised for HERA measure- 
ments of proton structure functions which are dominated by gamma exchange, therefore the 
parametrised PDFs are the valence distributions xuy and xd^, the gluon distribution xg, and 
the u-type and d-type sea quark distributions xJJ, xD, where xJJ = xu, xD = xd + xs. Us- 
ing a simple parametrisation form with the normalisation parameters constrained by the QCD 
sum-rules and relying on additional assumptions^ the number of free paramaters reduces to 10 
for the fits to HERA I data. However, the more data become available, the more constraining 
assumptions can be released, therefore the number of free parameters once HERA II data are 
added is increased, allowing for more freedom to the gluon and valence distributions, HER- 
APDFl.Sf. For example addition of the preliminary combined charm data, or jet data from 
HlEEl and ZEUSEEI allows for more fiexibility of the gluon distribution. The PDFs are then 
evolved using the DGLAP evolution equations ^ at NLO and NNLO in the MS scheme with 
the renormalisation and factorisation scales set to Q^. The QCD predictions for the structure 
functions are obtained by convoluting the PDFs with the calculable coefficient functions tak- 
ing into account mass effects for the heavy quarks based on the general mass variable flavour 
scheme'^. The PDF uncertainties at HERA are classified in three categories: experimental, 
model, and parametrisation. The consistency of the input data set and its small systematic 
uncertainties enables us to calculate the experimental uncertainties on the PDFs using the x^ 
tolerance Ax^ = 1. The model uncertainties are evaluated by varying the input assumptions, 
as performed in'^'. The parametrisation uncertainty is estimated as an envelope which is formed 
as a maximal deviation at each x value from the central fit. 

3 Results and Comparisons 

The inclusion of the precise high Q^ preliminary HERA II data in the QCD fits results in 
HERAPDF1.5. Figure[l]shows the comparison between HERAPDF1.5 and HERAPDFl.O which 
is based on HERA I data alone. The impact is noticeable especially in the high x region, where 
the valence contribution dominates. In addition, a meticulous study has been performed to 
estimate the uncertainties for the NNLO HERAPDF set, also shown in Figure [T| where the new 
NNLO set is compared to the HERAPDFl.O NNLO. To answer the question about implications 
of the new PDF sets on the Higgs exclusion limits from Tevatron, Figure [2] shows a more explicit 
comparison of the gluon distribution at high x between HERAPDF and MSTW08 sets at NNLO. 
The HERAPDF1.5 NNLO set yields a harder gluon at high x compared to the HERAPDFl.O 
NNLO set, which is prefered by the Tevatron jet data. The differences are mostly due to the 
use of a more flexible functional form made possible via the availability of a more precise data. 
Another important result is related to the determination of the strong coupling. New results 
from HERA based solely on the inclusive DIS measurements can not yet pin down the value of 
as as shown in Figure^ (right), where a scan in x^ ^^ function of strong coupling is shown. 
However, as soon as the jet data are included in the QCD flt, resulting in HERAPDF1.6, it 
can be seen that the strong coupling is well constrainted. In fact, the addition of the HERA 
jet data into the fit allows the strong coupling and gluon to be simultaneously constrained. 
Figure [3] shows the impact of the inclusion of the jet data on the gluon distribution when the fit 



is performed with the strong couphng as a free parameter. As a result, the HERA data prefer a 
rather larger value for as{Mz) = 0.1202 ± 0.0013(exp) ± 0.0007(mod) ± 0.0012(had)to:[}o3i(th). 

In addition, the inclusion of the preliminary HERA charm data provides constraints for the 
optimal value of the charm mass used in theory models. It has been observed that QCD fits 
without charm data have only a small sensitivity to the value of the charm mass. However, 
there is a strong preference for a particular rric once charm data is included. This study has 
been performed for various schemes, such as those used in the global fit analyses of MSTW08 
and CTEQ. The results conclude that each scheme describes the data well at the corresponding 
best value of the rric. It is interesting to observe that differences in the PDF sets correspond to 
differences in the charm mass used in different schemes'^. 

All HERAPDF sets derived solely from ep measurements are able to give a good prediction 
of the Z and W at Tevatron from the pp processes, and provide a competitive description to pp 
processes at the LHC as well. 



HI and ZEUS Combined PDF Fit 

Q- = 10 GeV^ 



HI and ZEUS HERA I+U PDF Fit 
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HERAPDF1.5 NNLO (prel.) 

^^1 exp. uncert. 

\^^ model uncert. 
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Figure 1: On the left hand side: the comparison between HERAPDF 1.0 (light colour) based on the HERA I data 
and HERAPDF1.5 (dark colour) based on the HERA I and H data, using the total uncertainty band at Q^ = 10 
GeV^ with gluon, sea (which are scaled by a factor of 0.05) and valence distributions. On the right hand side: the 
summary plot for HERAPDF1.5 at Q^ = 10 GeV^ at NNLO with the uncertainties including the experimental 
(red), model (yellow) and the PDF parametrisation (green), compared to the central fit for HERAPDFl.O at 

NNLO. 



4 Summary 



HERA provides accurate determinations of the proton structure and can predict related Stan- 
dard Model processes. Additional preliminary combined measurements from HERA II allow the 
high X region to be better constrained resulting in a more precise HERA PDF set, HERAPDF 1.5 
at NLO and NNLO in as with a detailed analysis of the uncertainties. The HERAPDF sets, 
which are based solely on ep data also describe the Tevatron data well and provide competitive 
predictions to the LHC processes. Inclusion of the HERA jet data allows for a precise deter- 
mination of the strong coupling, for which HERA data prefers a rather larger value of as- In 
addition, inclusion of the preliminary combined HERA charm data provides constraints on the 
optimal value for the charm mass used in theory models and it may account for some of the 
differences among global PDF fits. Therefore, HERAPDF provides a large variety of meticulous 
studies based on new measurements from HERA for accurate predictions at the LHC'^l. 
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Figure 2: On the left hand side: the gluon distribution from the HERAPDF sets (based on HERA I only - 

HERAPDFl.O, and based on HERA I+H, HERAPDF1.5) to MSTW08 at NNLO in the high x region. On the 

right hand side: the Ax^ distribution as a function of the value of as{Mz) in the PDF fits for HERAPDFl.Sf 

(dashed line) without jet data and HERAPDF1.6 (solid line) with jet data. 
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HI and ZEUS HERA I+II PDF Fit 
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HI and ZEUS HERA I+U PDF Fit with Jets 
Q' = 10 GeV^ 
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Figure 3: Summary plot for HERAPDF1.5f (left) and HERAPDF1.6 (right) as a function of x for Q^ = 10 GeV^ 
The strong coupling constant as{Mz) is a free parameter in both fits. The central values of the PDFs (solid 
lines) are shown together with the experimental, model and parametrisation uncertainties represented by the red, 

yellow and green shaded bands, respectively. 
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